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ABSTRACT 



Pluto's icy surface has changed colour and its atmosphere has swelled since its last 
closest approach to the Sun in 1989. The thin atmosphere is produced by evaporating 
ices, and so can also change rapidly, and in particular carbon monoxide should be 
present as an active thermostat. Here we report the discovery of gaseous CO via the 
1.3mm wavelength J=2-l rotational transition, and find that the line-centre signal 
is more than twice as bright as a tentative result obtained by Bockelee-Morvan et 
al. in 2000. Greater surface-ice evaporation over the last decade could explain this, or 
increased pressure could have caused the atmosphere to expand. The gas must be cold, 
with a narrow line-width consistent with temperatures around 50 K, as predicted for 
the very high atmosphere, and the line brightness implies that CO molecules extend 
up to w 3 Pluto radii above the surface. The upper atmosphere must have changed 
markedly over only a decade since the prior search, and more alterations could occur 
by the arrival of the New Horizons mission in 2015. 
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1 INTRODUCTION 

Pluto was the first body to be discovered in the Kuiper 
Belt of icy planetesimals orbiting in the outer Solar System. 
Although now only the second-largest body known in the 
belt (after Eris), Pluto is intensively studi ed, and is the pri- 
mary target of the New Horizons mission (|Stern &: Spencerl 
2003), with a spacecraft fiyby scheduled in 2015. In partic- 
ular, Pluto offers the rare opportunity to study a very cold 
planetary atmosphere which also evolves strongly with time, 
and may in fact be transient, settling and freezing out on 
the surface l|Owen et al.l 1 19931 ) as Pluto recedes from the 
Sun in its highly elliptical orbit. As Pluto's last perihelion 
was in 1989, subsequent contraction of the atmosphere was 
expected, but in fact when observed by occultation of stars 
the atmospheri c pressure and siz e had increased between 
1988 and 2002 (|Elliot et alj|2003ej ). Around the same time, 
the surface ice had become redder, and also brig hter around 
the no rthern pole as it came more into sunlight (|Buie et al.l 
l2010al lb[). Many of these results were surprising, such as 
changes in albedo over as little as a year, and darkening of 
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the southern regions, contrary to expectations if ices have 
been sublimating from here most recently. 

The atmosphere and surface must be strongly cou- 
led, since sublimat ion of ices supplies the atmosphere 
Schaller et al.l [20071 . e.g.), an d so as the com position of 
the surface changes seasonally (|Buie et al. 2010a, tj, the gas 
abundances shou ld also evolve. Atm ospheric temperatures 
of around 100 K ilElliot et alj|2 003b) are much warmer than 
the surface at only 40 K (|Trvka et all 119941 ). This requires 
a complex thermal balance controlled by the relative abun- 
dances of molecules that absorb radiation and re-emit in 
cooling lines, plus there is an overall influe nce of the leve l 
of sunlight, particularly in the ultra-violet (|Strobelll2008h . 
The atmosphere is also very fragile, with a s olar- driven slow 
wind predicted to be removing it into space jDelamerdl2009l . 
e.g.). 

Such a cold diffuse atmosphere is difficult to study 
by remote observation, and most conclusions are presently 
model-based. The trace gas methane is the only one t o 
have been confirmed spectroscopically (| Young et al.l I1997T ) . 
Surface-ice spectr oscopy shows frozen nitrogen dominates 
ijOwen et al.lll993T ). with metha ne and carbon mono xide also 
prese nt, plus possibly ethan e l|DeMeo et al.l [2O10l ) and ni- 
triles (jProtopapa et al1l2008l ). N2 must dominate the atmo- 



2 



J. S. Greaves et al. 



Radio velocity (USB) (km/s) 
-10455 




Radio velocity (LSB) (km/s) 

Figure 1. Spectrum towards Pluto in the J=2-l rotational transition of CO at 1.3mm wavelength. The atmospheric line from Pluto 
appears near zero velocity (in the dwarf planet's rest-frame) while other lines are residuals from background Galactic clouds varying 
from night to night. Summed over eight spectral channels, the Pluto line is detected at 6.5<r. 



sphere as it vaporiz es readily and is by a factor of 50 the 
most abundant ice (lOwen et al.l |l993 h . but sublimated CO 
is the strongest coolant l|Strobel|[200^ and so should act as 
the thermostat. Some evidence for CO includes t he brightest 
patch on Pluto's surface, suggested to be CO ice (|Buie et al.l 
l2010bl . and references therein), and an extinction layer seen 
in at mospheric occultation, interpreted as droplets of N2 or 
CO (jRannou fc Durry||2009r i. 



1.1 Previous searches for CO 

Prior searches for gaseous CO have been made with ground- 
based telescopes. Infrared absorption spectroscopy has re- 
cently constrained C O:N2 to be <0.5 % near the surface 
(iLellouch et alJl201Ch an order of magnitude below an ear- 
lier limit ( Young et al.ll200lT ). There is also a long history of 
millimetre-wavelength searches, looking for ro tational tran- 
sitions of CO from the hig her atmosphere (|Barnes! 1 19951 ; 
iBockelee-Morvan et aljfeooil ). As the atmosphere is warmer 
than the surface, and extends beyond the planetary disc, 
emission lines are predicted. The major difficulty is that 
Pluto is small and distant, subtending less than a percent 
of the beam solid angle, even for telescopes of 10-30m class. 
An upper limit was first obtained using the Haystack 37m 
in 1993, for the CO J=l-0 tr ansition rotational line at 2.6 
mm wavelength |Barnei[l993). This limit was improved by 
the equiv alent of eight-fold in 2000, w ith the IRAM 30m 
telescope (|Bockelee-Morvan et alj |2001). and this group also 
searched for the J=2-l transition at 1.3 mm with better sen- 
sitivity. A tentative positive signal was seen, with the bright- 
est spectral channel observed at 28 mK in main-beam bright- 



ness temperature, versus la noise of 10 mK per 0.1 km/s 
spectral channel. However, the observation was severely af- 
fected by a strong line from a background Galactic cloud, 
coming as close as 0.5 km/s to potential Pluto emission 
|Bockelee-Morvan et ai1l200ll ). 

The 2000 data suffered from an unfortunate coincidence 
in sky position with an uncatalogued Galactic source, but 
Pluto's orbit as seen from the Earth has brought it closer 
to the Galactic Plane over the last couple of years. Here 
we present results from the 15m James Clerk Maxwell Tele- 
scope, making a new search for the J=2-l line, at epochs 
chosen to minimise Galactic contamination behind Pluto. 
The search was motivated by time availability at an excellent 
millimetre-site, and the opportunity to examine the evolu- 
tion of Pluto's atmosphere i n the exciting run-up p eriod to 
the arrival of New Horizons ()Stern fc Spencerll2003l ). We re- 
port the confirmation of atmospheric CO (after two decades 
of searches) , making it only the second gas-pha se species to 
be detected, after methane ()Young et al J 119971 ) . 



2 OBSERVATIONS AND REDUCTION 

We used the JCMT at 4000 m altitude on Mauna Kea, 
Hawaii, observing over 3 nights in August 2009 and 8 nights 
in April to May 2010. Conditions were good to average for 
the site, with zenith sky opacity at 1.3 mm typ ically around 



0.1. T he setup used the ACSIS spectrometer (Buckl e et al 
2009) together with the A3 receiver (|Cunningham etal 



1992) (the instrument described was subsequently rebuilt by 
the Herzberg Institute of Astrophysics with a 1.3 mm-band 
mixer). Total on-source integration times were 12 and 19 
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Figur e 2. Left: expanded version of Figure 1. Overplotted in red is the spectrum obtained in April 2000 by iBockelee-Morvan et al,l 
(2001) using the IRAM 30m, but re-scaled for equivalent signal with the 15m JCMT; the brightest channel near zero velocity is then 
7± 2.5 mK over 0.1 km/s. Right: comparison of the average spectrum and the 2010 and 2009 data (top to bottom, offset for clarity). 
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Figure 3. Example spectrum if all 11 nights of data are weighted 
and co-added but no baselines are removed. The red channels were 
used to fit an example linear baseline, and this would extrapolate 
across the Pluto velocities as shown by the red straight line. 



hours in 2009 and 2010 respectively, with individual observa- 
tions lasting approximately 20 minutes. The observing mode 
was frequency- switching, with a throw of 35.5 km/s within 
the passband. The telescope tracked Pluto's sky position 
but not its velocity, so spectra were shifted in software to 
the dwarf planet's reference frame. The spectrometer chan- 
nel spacing was set to 0.04 km/s and there are systematics of 
around ±0.025 km/s for velocity-drifts of Pluto within each 
observation that have not been corrected. Data are shown on 
the main-beam antenna temperature brightness scale, with 
a beam efficiency of 0.75 at 230.538 GHz frequency. 

For each night of data, the shifted spectra were co- 
added, and then the appropriate positive and negative 
frequency-switched features were averaged. The region of the 
net positive line was blanked and replaced with a linear in- 
terpolation of surrounding channels. A running average over 
12 channels was then generated from this and subtracted 
from the pre-blanked spectrum as a baseline. The 11 nightly 
spectra made in this way were finally combined, weighted 
by 1/noise 2 factors to take into account their different con- 
ditions and observing durations. Rejected data include two 
earlier nights using a smaller frequency-switch, and two later 



nights when Pluto had approached within 1.5° of the Galac- 
tic Plane. These spectra show ripples and complex baselines, 
degrading the final line profile. The weight of these rejected 
data represent only 15 % of the total observations made. 

The observations were planned to maximise the height 
of Pluto from the Galactic Plane, which varies as seen from 
the Earth, to minimise contamination problems. Pluto's 
Galactic latitude in 2009 was +2.2°, and in 2010 was -2.0° 
to -1.7°, with longitudes of ~ +12°. Using two periods well 
spaced around the year (Aug and Apr/May) also caused 
large changes in Pluto's velocity relative to the telescope, 
from +26 km/s in 2009 to -27 to -19 km/s in 2010. Thus 
any narrow line seen persistently at Pluto's velocity in the 
co-added data can only be from the planetary atmosphere, 
as any Galactic background emission is different every night. 
Galactic lines were seen around Pluto's velocity in the spring 
data, but were much broader than the Pluto line. 



3 RESULTS 

Figure 1 shows the co-added Pluto spectrum from 2009/10, 
in the dwarf planet's rest-frame. A clear line is seen at zero 
velocity, with 6.5a confidence in the integrated antenna tem- 
perature. There are no other significant narrow features; the 
next brightest (at -5 km/s), is at only 2.7ct over three spec- 
tral channels, and was traced to a signal from a pair only 
of the 11 observing nights. The broad features are resid- 
uals from the running-average subtraction process, arising 
from broad lines of Galactic clouds from different nights. 
A detail spectrum (Figure 2) shows the very flat baseline 
achieved near Pluto's velocity, while a test made by blank- 
ing a broader region (20 channels instead of 12) showed a 
change in integrated brightness of only around 10 %. This 
argues for a robust signal measurement of the narrow line 
seen, and also implies that broad line wings are faint. 

The CO line from Pluto was independently detected in 
both the 2009 and 2010 epochs, at levels of 3.9er and 5.4a re- 
spectively, and with the same integrated brightness within 
the errors (Figure 2, right panel). The antenna tempera- 
tures are calibrated to about 20 %, based on rms scatter of 
the line-peak signal of spectral standards taken every night 
on the nearby source OH17.7-2. Also, inclusion of the two 
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nights of early data taken with a smaller frequency-switch 
was found to reduce the final integrated line intensity by 
« 20 %. The brightness measurements could be affected by 
broad structure in the baseline, and Figure 3 demonstrates 
the spectrum obtained if the data are co-added without 
baselining. The curvature due to Galactic residuals makes 
it difficult to establish the nature of faint line wings. 

The integrated antenna temperature of CO J=2-l over 
8 spectral channels (0.32 km/s) is 6.3±1.0 mK km/s on a 
main-beam brightness temperature scale. The mean T m j, is 
20 mK, and the la noise p er 0.04 km/s spectral channe l 
is 8.6 mK (twice as deep as IBockelee-Morvan et al.l (l200ll ) 
for a common spectral resolution). The full- width at half- 
maximum intensity is approximately 0.28 km/s. The ve- 
locity centroid is ~ +0.07 ± 0.055 km/s, including errors 
from comparing adjacent half-line intensities of 0.03 km/s 
and systematics of the 0.04 km/s channel spacing and 
±0.025 km/s velocity drifts within scans. Hence, the centroid 
velocity may be marginally positively offset from Pluto. 

Figure 2 also shows the CO J=2-l sp ectrum obtained 
in 2000 by IBockelee-Morvan etTafl (|200ll V Their T mb val- 
ues have been scaled downwards by a factor of four, to ac- 
count for the greater beam-dilution for a compact source 
with the JCMT 15m versus the IRAM 30m telescope. On 
this 'JCMT-equivalent' scale, the 2000 spectrum has a 3a 
upper limit of 7.5 mK in each 0.10 km/s channel, and the 
brightest channel near zero velocity is 7 mK. Weak positive 
signals were seen in line wings, and the integrated signal 
over ±1.3 km/s was 5.9 mK km/s, at signal-to-noise ratio 
of 4.5. The integrated intensity appears similar in 2000 and 
2009/10, but the earlier spectrum has a greater contribution 
from possible line wings (Figure 2). A Voigt-profile fit to the 
new spectrum suggests broad wings now contribute < 35 %. 



3.1 Features of the CO spectrum 

Figure 2 shows that the line-centre emission in 2009/2010 
significantly exceeds tha t in th e observation from April 2000 
jBockelee-Morvan et all l200ll ). In the brightest 0.1 km/s 
channel, the signal was 7 ± 2.5 mK, compared to 27 ± 5 mK 
now when averaged over the same region from the JCMT 
spectrum. Thus the ratio is at least 22/10.5 from the la 
bounds, implying the line-centre intensity has more than 
doubled over a decade. 

The core of the line now seen is narrower than 
previously p r edicte d in models |Stansberrv et all 1 19941 ; 
IStrobel et alj fl996l). Model line profiles as plotted by 
IBockelee-Morvan et alj (|200ll ) are over twice as wide as 
observed, with full-width half-maxima of approximately 
0.7 km/s. These models had atmospheric temperatures of 
80-106 K, but our measured line width of 0.28 km/s implies 
that thermalised CO can not be warmer than ~ 50 K. This 
estimate includes only thermal broadening, and so the gas 
could be e ven cooler if t here is another source of velocity 
dispersion. IPerson et al.l (120081 ) find evidence for large-scale 
atmospheric waves, but with horizontal wind speeds of under 
3 m/s, wh ile upwards sublim ation flows should be at only a 
few cm/s (|Toigo et alj|20ich . The latter would also be in the 
contrary sense to the tentative red-shift noted above, since 
upwelling would produce a blue-shift from the near-side of 
the atmosphere that is preferentially observed. 



4 DISCUSSION 

The marked increase in CO line-peak brightness could be 
due to recent sublimation of surface ice - for example, the 
largest bright patch see n on Pluto's surfac e is attributed to 
carbon monoxide frost (|Buie et al.ll201o"bh . While this is a 
long-term feature in albedo maps, it could be a rich source of 
sublimating CO. Alternatively, a large expansi on in isobar 
height s was detected between 1988 and 2002 (|Elliot et al.l 
l2003al ). and such increases in atmospheric pressure could 
raise the effective emitting area for trace molecules. How- 
ever, these occultati on data only detect the atmosphere up 
to around 135 km l|Young et alj [20081. e.g.), so effects at 
higher altitudes are uncertain. Regarding timescales, marked 
surfa ce changes have o ccurred over intervals as short as a 
year (|Buie et alj|2010bf ). while vertical sublimation speeds 
l|Toigo et alj|2010l ) imply a molecule could rise thousands of 
kilometres within a decade. Hence it is reasonable that tem- 
poral changes in gas densities and/or abundances as surface 
ices sublimate could result in the CO lines being quite dif- 
ferent over observations a decade apart (Figure 2). 

Finding cold CO points to gas at very large heights - for 
example temperatures as low as 50 K have been predicte d 
for about 4500 km above the centre of Pluto l|Strobelll2008h . 
This is under solar minimum conditions, as actually present 
in 2009/10. If the top of the CO layer is at 4500 km, then 
the atmosphere would subtend a diameter of 0.4 arcsec, as 
it was seen in 2009/10 with Pluto at 31 AU from the Earth. 
From this, in the simplest approximation of a filled opaque 
shell of CO molecules thermalised at 50 K, the beam-filling 
factor would be 3.6x 10 -4 , giving an expected line brightness 
of 18 mK. This is in good agreement with our mean T m t of 
20 mK. Lower atmospheric layers should be warmer and so 
add line brightness, but as an example, ga s twice as ho t 
(100 K) is predicted to lie below 2000 km (|Strobel|[2008l ). 
and thus be five times more beam-diluted. 



4.1 Comparison to models 

Exact calculations for the CO line profile require radiative 
transfer calculations and an atmospheric model, beyond the 
scope of this Letter. However, some estimates can be made 
for the heights of the emitting layers. The density profil e 
should decline approximately as a power-law i|Strobelll2008h . 
if CO has a constant mixing ratio with N2. Then the height 
below which the CO emission becomes opaque can be es- 
timated, or more strictly, the height of the tangent where 
the column of CO molecules excited to the J=2 level pro- 
duces J rdv ~ 1 x Aw, for opacity r of unity and Av equal 
to the line FWHM of 0.28 km/s. In an isothermal 50 K 
approximation, this height is 2750 km for a power-l aw den- 
sity p rofile approximating that at solar maximum (|Strobell 
2008), and 2250 km for a solar minimum profile. The respec- 
tive heights which contribute opacity of 0.01 are about 4600 
and 3200 km. The se values are for CO:N 2 of 0.005, follow- 
ing IStrobell (l2008h and at an upper limit observed in July 
2009 l|Lellouch et al.ll2010l) . If contributions to the CO signal 
come from as high up as 4500 km, as estimated above from 
beam-filling arguments, then it appears that atmospheric 
densities must be at the high end of the models, or upper- 
level CO must be more abundant than assumed. Also, in the 
models high densities are associated with warmer temper- 
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atures than infered here - this anomaly may be related to 
heating from CH4 molecules absorbing solar U V, noting that 



recent estimates of the methane mixing ratio JZalucha et al 

mare up to ~ 5 times lower than assumed by IStrobe] 

The CO upper boundary estimated to be at about 
4500 km above the cen tre of Pluto correspon ds to ~ 4 radii, 
for iipiuto ~ 1153 km ijTholen fe Buidll997h . This is a rea- 
sonable upper bound, given that above ~ 3-5 Rpi uto, the 
gas s hould become so diffuse that is is collisionless (|Strobell 
I2008T 1. and thus the CO 3=2 level will not be populated. 
Also, interaction with the solar wind will cause ionised 
gas (and any swept-up neutrals) to flow away into space , 
above sa 5 R p iuto i n the simulations of iDelamerd (2009); 
lHarnett et al] (|2005l ) . The CO layer estimated to extend out 
to ~ 4 Rpiuto is t hus inside this flow zone and can be ex- 
pected to be stable. IStrob el (2008) models a slower hydrody- 
namic outflow inside this region, but predicted speeds only 
exceed our 0.04 km/s velocity resolution above 4500 km, so 
effects on the spectrum would be negligible. It is however 
possible that the marginal CO line red-shift, if real, could 
indicate a flow forming into a comet-like tail directed away 
from the Sun. 

CO gas should be chemically stable even at large alti- 
tudes, and at low densities, cloud models imply that it will 
stay in the gas phase. CO ice could only form if seed par- 
ticles existed - akin to noctilucent clouds on Earth, where 
highly supersaturated water can condense out on rare solid 
particles, such as meteoritic dust raining in from space. Such 
interp lanetary dust is how e ver sparse out in the K uiper 
Belt (|Landgraf et alj l2002h . IR annou fe Durrvl (|2009l ) pro- 
posed cloud decks lower down in Pluto's atmosphere, with 
droplets of CO or N 2 that could explain an extinction layer 
ijElhot et all2003al ). However, heterogeneous nucleation (i.e. 
chemical reactions on a surface leading to increase of the 
droplet volume) is only efficient inside a limited tempera- 
ture interval, and it yields no droplets above heigh ts of ~30 
km in these Pluto models (|Rannou~fe Durry||2009h . 

Since surface changes dBuie et al 1 l2010alibh and the 
gaseous methane abundance (|Zalucha et al.ll201ll ) are seen 
to be temporally variable, time-dependent models will be 
needed to balance the effects of heating, cooling, sublima- 
tion and solar irradiation. Further CO line monitoring is un- 
derway, and N2 gas can be detected for th e first time with 
ultraviolet spectroscopy (IStark et all 120071') b y the ALICE 



b y tne 

spectrograph on New Horizons (jStern et all 120081 ). Hence 
atmospheric models can be tested from data taken before 
and during the 2015 fly by. 



5 CONCLUSIONS 

Carbon monoxide is confirmed in the atmosphere of Pluto, 
and gives the first detection of high-altitude gas. There have 
been large temporal changes in the CO line profile over the 
last decade, in parallel to major alterations of the surface 
ices. Pluto's rapidly changing and fragile atmosphere could 
provide an interesting test-bed for models of global climates. 
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